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Abstract 



We present new analytic and numerical results for self-gravitating SU(2)- 
Higgs magnetic monopoles approaching the black hole threshold. Our investi- 
gation extends to large Higgs self-coupling, A, a regime heretofore unexplored. 
When A is small, the critical solution where a horizon first appears is extremal 
Reissner-Nordstrom outside the horizon but has a nonsingular interior. When 
A is large, the critical solution is an extremal black hole with non-Abelian hair 
and a mass less than the extremal Reissner-Nordstrom value. The transition 
between these two regimes is reminiscent of a first-order phase transition. We 
analyze in detail the approach to these critical solutions as the Higgs expec- 
tation value is varied, and compare this analysis with the numerical results. 
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I. INTRODUCTION 



The physics of a nonsingular spacetime is qualitatively distinct from a that of a space- 
time exhibiting a black hole. However, families of spacetimes exist that may be viewed as 
interpolating between the two. These spacetimes are nonsingular and have no horizons; 
nevertheless, they have a region whose metric can be made arbitrarily close to that of the 
exterior region of a black hole. In the limiting case, the inner boundary of this region takes 
on the characteristics of an extremal horizon, even though no curvature singularity develops 
in the interior. By studying such solutions, one may gain further insight into the properties 
of black holes. 

One approach to the construction of such "almost black holes" begins with a sponta- 
neously broken Yang-Mills theory that has magnetic monopole solutions. For small Higgs 
expectation values, these monopole solutions persist when gravitational effects are included. 
However, as the Higgs expectation value is increased toward a critical value on the order of 
the Planck mass, the monopole solutions begin to approximate black holes with finite mass 
and nonzero horizon radius. 

Such gravitating monopoles, as well as the related magnetically charged black holes with 
hair, have been studied previously JD-EI]. For a review article and recent work in related 
subjects, see ||. In this paper we investigate these solutions in more detail, concentrating on 
some aspects that were not previously noted. Most notably, we find that two distinct types 
of behavior, with qualitatively different extremal black hole limits, can occur. Our focus 
here is on the detailed properties of these solutions and their behavior as they approach the 
black hole limit. We describe elsewhere how these solutions can be used as approximate 
black holes that provide insight into the transition from a nonsingular spacetime to one with 
horizons. 

We work in the context of an SU(2) gauge theory that is spontaneously broken to U(l) 
by a triplet Higgs field with vacuum expectation value (0) = v. The elementary particle 
spectrum of this theory contains a pair of vector mesons that carry U(l) electric charge 
and have mass ev (where e is the gauge coupling), as well as a massless photon and a 
massive electrically neutral Higgs particle. In flat space the classical field equations have 
a spherically symmetric monopole solution with U(l) magnetic charge 4n/e and a mass of 
order Airv/e. This has a central "core" region, of radius -R cor e ~ ^-l^w, in which there are 
nontrivial matter fields. Beyond this radius is a "Coulomb" region in which the massive 
fields fall off exponentially fast, leaving only a long-range Coulomb magnetic field. 

In studying the behavior of these solutions in the presence of gravity, we assume spherical 
symmetry, and so can write the metric in the form 



ds 2 = B(r)dt 2 - A(r)dr 2 - r 2 (d6 2 + sin 2 9d<j) 2 ) . 




It is often convenient to rewrite A in terms of a mass function m(r) defined by 
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with m(oo) = M. For a configuration to be nonsingular at the origin, A(0) = 1 and 
m(0) = 0. A horizon occurs when 1/A has a zero or, equivalently, when m{r)/r = 1/2G. 

A benchmark with which to compare our results is provided by the Reissner-Nordstrom 
metric, with 



B(r) = - = 1 - — + ^ (13) 
A r Anr 2 

If M > \J Q 2 1 4:7rG , this describes a black hole solution with a charge (either magnetic or 
electric) Q and an outer horizon determined by the larger of the two zeroes of 1/A. If instead 
M < d l Q 2 1 '4nG , there is no horizon separating the curvature singularity at r = from 
the asymptotic regions and a naked singularity results. The boundary between these two 



regimes, M = y Q 2 / AirG , gives the extremal Reissner-Nordstrom black hole. For Q = Air/e, 
the case with which we will be concerned in this paper, the extremal black hole has a horizon 
radius 



ro = \/-5- ( L4 ) 



e 2 



and a mass 

M =g. (1.5) 
When v <C Mp\ the gravitational effects on the monopole are relatively small. The metric 



at large distances approaches the Reissner-Nordstrom form, with \JQ 2 /AnG M. There is 
no singularity, because the actual metric deviates from the Reissner-Nordstrom form when 
r -Rcore- One finds that 1/A = 1 at the origin, decreases to a minimum at a radius of 
order R COTe , and then increases monotonically. As v is increased, the core shrinks and the 
minimum of 1/A moves inward and becomes deeper. In most cases, this continues until 1/A 
develops a double zero, corresponding to an extremal horizon, when v = v CT . A slightly 
different behavior is found for very small Higgs self-coupling [^]. In this case, the solution 
varies continuously as v is increased up to a value f max . Although the minimum value of 
1/A is still nonzero at this point, static solutions do not exist for higher values of v. Instead, 
these solutions join smoothly on to a second branch of solutions for which (1/A) m i n decreases 
as v is decreased from w max to a critical value v CI , where the extremal horizon develops. 

For small values of the Higgs mass, the extremal horizon of the critical solution occurs at 
the Reissner-Nordstrom radius tq, in the exterior Coulomb region of the monopole solution. 
The matter fields take on their vacuum value everywhere outside the horizon and the exterior 
metric is exactly that of an extremal Reissner-Nordstrom black hole. When the Higgs to 
vector mass ratio is greater than about 12, a regime not explored in depth in previous studies, 
the behavior is rather different. The horizon in the monopole core at a radius r* < r that 
decreases with increasing Higgs mass. In this case there are nontrivial matter fields, or 
"hair", outside the horizon. The transition between the two regimes is not smooth, but 
instead is reminiscent of a first order transition. 
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The remainder of this paper is organized as follows. In Sec. II we outline the general 
formalism and define our conventions. In Sec. Ill we use numerical methods to obtain 
monopole solutions to the field equations. We describe in detail their behavior as v is 
increased towards its critical value. The critical solutions that are the limits of these families 
of monopole solutions are characterized by the presence of extremal horizons. In Sec. IV, 
we use analytic methods to study the properties of these critical solutions, focusing on the 
behavior near the horizon. We show that the problem of finding a solution with an extremal 
horizon can be formulated as a pair of boundary value problems, one for the region < r < 
and one for r* < r < oo, that must be solved simultaneously. The conditions for a solution 
to these places strong constraints on the behavior of the fields near the horizon. These 
constraints allow only two types of behavior, one associated with a core region horizon, the 
other with a Coulomb region horizon, which we examine in detail in Sec. V. In Sec. VI we 
compare these analytic predictions with our numerical results. Section VII contains some 
concluding remarks. The Appendix contains details of the numerical investigation of the 
transition region between the two types of critical solutions. 



II. GENERAL FORMALISM 



We consider an SU(2) gauge theory with a triplet Higgs field a whose self-interactions 
are governed by the scalar field potential 



(2.1) 



In flat spacetime, this theory has nonsingular monopole solutions, with magnetic charge 
Aw /e, that are described by the spherically symmetric ansatz 



cj) a = vr a h(r) 
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(2.2) 

(2.3) 
(2.4) 



While it is clear that vh(r) is the magnitude of the Higgs field, the meaning of u(r) is some- 
what obscured by this "radial gauge" ansatz. By applying a singular gauge transformation 
that makes the direction of the Higgs field uniform, one finds that u(r)/er is equal to the 
magnitude of the massive vector field, and so it, like 1 — h(r), should be expected to vanish 
exponentially fast outside the monopole core. 

The generalization of this ansatz to curved spacetime is straightforward. Since we are 
considering only static spherically symmetric solutions, we can take the metric to be of the 



form of Eq. (1.1). The matter part of the action is 



^matter = -47T J dtdrr 1 \f~AB 



K(u,h) 
A 



+ U{u, h) 



(2.5) 
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where 

K 



2 2 / „ x 2 



1 /duV v 2 dh' 
e 2 r 2 \dr J 2 \ dr ) 
l u 2 - l) 2 v 2 u 2 h 2 Xv 4 

2e 2 r 4 + r 2 

One may view U(u, h) as a position-dependent potential. It has several stationary points, 
which we enumerate here for later reference: 



1. u = 0, h = ±1. This is a local minimum of C/ for r > 1/ev. 

2. u = u(r), h = h(r), or u = —u(r), h = —h(r), where 



/ A(l - e 2 v 2 r 2 ) 

M= V X-e 2 (2 ' 6) 



^J A 7 2 _ /( e r )2 - (2.7) 



These are real only for r lying between 1/ew and 1/a/A In that range it is a minimum 
of U if A/e 2 > 1, but a saddle point if X/e 2 < 1. If A/e 2 = 1, this solution is replaced 
by a degenerate set of local minima, with h 2 + u 2 = 1, that exist only for evr = 1. 

3. u — h — 0. This is never a local minimum of U. 

A. u = ±1, h = 0. This is a local minimum of U for r < l/\/Av. 

We will see that these stationary points are key to understanding the local existence of 
extremal horizons in monopole systems. 

For static solutions, the matter fields obey the equations 



1 d fVABdu\ e 2 r 2 dU u(u 2 - 1) 2 2 l2 

1 1 - 7 +eVw/z 2 (2.8) 



VAB dr \ A dr ) 2 du r 2 
1 d fr 2 VABdh\ 1 917 2/m 2 



r 2 \/AB dr \ A dr J v 2 dh r 2 
These must be supplemented by the two gravitational field equations 



+ 2Xv 2 h(h 2 - 1) . (2.9) 



1 dVJB =8,GrK (2.10) 



VAB dr 
d (\ 



l-8nGr 2 u) - ^-{l + 87iGr 2 K) . (2.11) 



r — , , , . , 

dr \AJ v > A 



Note that, up to a rescaling of distances, the solutions of these equations depend only on 
the dimensionless parameters^ a = 8itGv 2 and b = X/e 2 = {the /2mw) 2 ■ 



1 These parameters are related to those in Ref. Q by [i = a, in Ref. Q by a = 2b and (3 = a, and 
in Ref. [|| by a = yf^j2 and /3 = 2\fb. 
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Integration of Eq. (|2.10 ) gives B(r) in terms of the remaining functions. With the 
boundary condition A(oo)B(oo) = 1, corresponding to the conventional normalization of t, 
we have 

1 



Br) 



A(r] 



cxp 



16ttG / dr'r'K 



(2.12) 



Using this result to eliminate B(r) from the remaining field equations leaves one first-order 
and two second-order equations to be solved. A solution of these is determined by five 
boundary conditions. Requiring that the fields be nonsingular at the origin gives three of 
these, w(0) = 1, h(0) = 0, and A(0) = 1. Two more, w(oo) = and h(oo) = 1, follow from 
the finiteness of the energy. Additional boundary conditions arise when a horizon is present. 
While these are not relevant for our numerical solutions, which are all regular monopoles, 
they play an important role in our analysis of the extremal black hole configurations that 
are the limiting points of sequences of monopole solutions. 



III. MONOPOLE SOLUTIONS 



A. The method 



In our search for regular solutions and their approach to criticality, let us identify an en- 
ergy functional on the space of static configurations. Following van Nieuwenhuizen, Wilkin- 
son, and Perry ||, one can write down an action which is a functional of just u{r) and h(r) by 
not only eliminating V AB using Eq. ( |2.10|) , but also eliminating A by integrating Eq. ([2.1 1|) 
subject to the boundary condition A(0) = 1. The complete action (i.e., gravitational plus 
matter) then becomes 



S 



dt E 



j r d 



E = 4vr / dr { ——^-VAB 
/o \ SttG dr 



+ r 2 VAB 



1 ■ 

A. 
K(u,h) 
A 



(3.1) 



+ U(u, h) 



with v AB and A understood to be implicit functionals of u(r) and h(r). Configurations 
that extremize this action are solutions to Eqs. ( ET^ -( [2~7lT|) , with the energy E being equal 
to the mass M = m(oo) |§. This energy functional is not bounded from below, even when 
satisfying the appropriate boundary conditions, unless the 1/A{r) that corresponds to a given 
matter field configuration {u(r),h(r)} is everywhere greater than zero. This restriction is 
not serious for our purposes, since we are only interested in static regular monopole solutions 
and their approach to criticality; all such solutions satisfy the condition that 1/A(r) > 
everywhere. 

Thus, we find regular solutions to the field equations by seeking extrema of this energy 
functional that satisfy the boundary conditions at r = and r = oo. This is equivalent to 
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solving 5E/5u = 5E/5h = 0. We do this numerically by solving the alternative system of 
equations 

d 2 u ^du 5E 

dt 2 dt 5u 
d 2 h dh_5E 

dt 2+L dt~5h [6 - Z) 

for an appropriate damping factor F and taking the final steady state solution as the result. 
This approach is analogous to having a massive particle roll on a manifold determined by 
E while the particle's motion is viscously damped, eventually coming to rest at some local 
minimum of the energy. 

The discretization used to numerically implement this algorithm places a limit on how 
close we can approach the critical solutions in which the extremal horizon has actually 
formed. The errors in the fields are proportional to the square of the spatial step size. As 
a result is that we can only obtain solutions in which the minimum value of 1/A is at least 
0[(Ar) 2 }. 

One should contrast our relaxation method with the algorithms used in previous analyses 
]1|-[3|] that employ shooting from the origin. In these, a choice is made for the values of the 
fields and their first spatial derivatives at the origin, and the equations of motion are then 
used to integrate out towards infinity. The initial choice is then adjusted to ensure that the 
boundary conditions at infinity are satisfied. With this approach, difficulties appear with 
large b because of the extreme sensitivity at the origin to small perturbations. No such 
problem exists with relaxation, and as a result we have been able to obtain solutions for 
large b. However, although solutions that are unstable under perturbations may be found 
by shooting, such solutions cannot be obtained by our relaxation method. 

In particular, previous work || has shown that for b & 0.1 monopole solutions exist for 
a range of values of a that are greater than the value a cr of the critical solution. Moreover, 
for each value of a in this range there are two solutions, with the one having the smaller 
value of (1/A) min being unstable Hence, for this range of b our methods cannot find 

the critical limit of the regular monopole solutions. 



B. Approach to criticality: low —b case 



We investigate the approach to criticality by studying the nature of the monopole so- 
lutions as a = 8ttGv 2 is increased with the Higgs self-coupling held fixed. Two distinct 
behaviors are seen, depending on the size of b = X/e 2 . We will describe the low-6 behavior 
in detail in this subsection, and the high-6 behavior in the next. Figure [IJ illustrates the 
behavior for b = 1.0, a typical low-fe case. As a increases and gravitational effects become 
stronger, 1/A(r) begins to dip down until at a = a CT it develops a double zero, corresponding 
to an extremal horizon, at r = r , the horizon radius of the extremal Reissner-Nordstrom 
solution. At the same time, the matter fields u(r) and h(r) are pulled inward. At a = a cr , 
the variation of these fields occurs entirely within the horizon. Only the Abelian Coulomb 
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FIG. 1. Monopole solutions for b = 1.0 and various values of a. The progression from 
dot-dashed line, to dashed line, to dotted line, to solid line corresponds to a = 1.0,2.0,2.3, and 
2.37. The first three panels depict the metric functions (a) 1/A(r), (b) B(r), (c) (AB) 1 / 2 ^). (d) 
Matter fields variables; u(r) begins at unity at the origin and asymptotes to zero as r — > oo. h(r) 
is zero at the origin and asymptotes to unity as r — > oo. 
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versus (l/A) m - m for b = 16.0 and various values of a near o, 



magnetic field survives outside the horizon, while the metric for r > r$ is precisely that of 
the extremal Reissner-Nordstrom black hole. 

All of this reproduces results found in earlier work [!]-§]. Some features that were not 
previously stressed are revealed when we examine B(r). Coming in from large r, B decreases 
with 1/A until the latter reaches its minimum. B then continues to decrease, although at a 
much smaller rate. For very small h this decrease continues all the way in to r = 0, while 
for somewhat larger values of b there is a minimum in B at a finite r < r . The situation 
for the critical solution is somewhat ambiguous. If we adopt the conventional normalization 
B(oo) = 1, then B(r) vanishes identically inside the horizon. If instead we set B(0) = 1, 
then B is finite and varying inside the horizon and infinite outside the horizon, with a 
minimum either at r = or at some finite radius, depending on the value of b. In neither 
case is the minimum a zero of B. 

Closely related to this is the behavior of V ' AB. When a is small, this is very nearly 
constant, with a value close to unity. As a is increased, V AB develops a step-like behavior 
until, at criticality, it is precisely a step function centered at the horizon. This is in sharp 
contrast with the Schwarzschild and Reissner-Nordstrom solutions, where V AB = 1 every- 
where. To show how this behavior becomes more pronounced as the solution approaches 
criticality, in Figure |]we plot the value of V AB at the origin as a function of (1/A) min for 
b = 16. Note that there is a power law relationship between the two. Similar curves, but 
with different powers, are found throughout the low-6 regime. 

As a final illustration of the behavior of the low-6 solutions near criticality, in Fig. |3] we 
plot u for several near-critical solutions. However, rather than using the variable r, we have 
plotted XL cLS cL function of the proper length 
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1 

FIG. 3. A^ 1 (curves with minima) and u (monotonic curves) as functions of the proper length 
/ for b = 16.0. The progression from dashed line, to dotted line, to solid line corresponds to 
a = 1.529, 1.53 and 1.531 with a cr 1.532. 

l(r) = £ dr ^A(fj . (3.3) 

We see that, for sufficiently large I, u(l) is close to a decaying exponential, and that this 
behavior does not significantly change as one crosses the minimum of 1/A; similar behavior 
is seen with 1 — h. 

C. Approach to criticality: high —b case 

The scenario when b is large differs significantly from that for small b; we illustrate this 
in Fig. |] for the case of b = 100. So long as a is not too near its critical value, the qualitative 
evolution of the metric functions and field variables is similar to that in the previous case: 
1/A(r) dips down in a region r « r , the value of B(r) decreases in a region near the origin, 
and the fields u, h are drawn into the core of the system. 

Near criticality, a different behavior emerges. The evolution of 1/A(r) as the solution 
approaches criticality is depicted in Fig. |]a and, in more detail, in Fig. |^. Although initially 
similar to that the for the low- 6 case, the decrease of the minimum at r « r ceases before 
it actually reaches zero. As this occurs, a second minimum at a radius r < r rapidly 
drops down and forms a double zero, corresponding to an extremal horizon, at r = r* < To- 
Figure (and, in more detail, Fig. |6]) shows the function B(r) as a approaches its critical 
value. In contrast to the small-6 case, in the critical limit the minimum of B is a zero located 
at the horizon, r = r*. 

The evolution of V AB is shown in Fig. (|c. The behavior is similar to that for the small- 
b case until the outer minimum of 1/A stops decreasing and the inner minimum begins to 
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FIG. 4. Monopole solutions for b = 100 and various values of a. The progression from 
dot-dashed line, to dashed line, to dotted line, to solid line corresponds to a = 1.0, 1.1, 1.2, and 
1.29344. The first three panels depict the metric functions (a) 1/A(r), (b) B(r), (c) (AB) l l 2 (r). 
(d) Matter fields variables; u{r) begins at unity at the origin and asymptotes to zero as r — > oo. 
h(r) is zero at the origin and asymptotes to unity as r — > oo. 
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FIG. 5. Details of the metric function 1/A(r) near criticality at b = 100. The progression from 
dot-dashed line, to dashed line, to dotted line, to solid line corresponds to a = 1.27, 1.29, 1.293, 
and 1.29344. 
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FIG. 6. Details of the metric function B{r) near criticality at b = 100. The progression for 
different a's is the same as in Figure 5. 



11 



appear. Until this point, there is a power-law relationship between V AB\ r= o and (1/A) min 
similar to that for the small-6 case, with most of the variation in yAB occurring at r « r$. 
Once the inner minimum begins to appear, the decrease in V AB\ r= o slows down, so that 
even for the critical solution \J AB\ r=0 ^ 0. 

Figure f|d depicts the matter fields u{r) and h(r) for a series of values of a at fixed 
large b. Their behavior is analogous to that for small b except that the fields are never 
completely drawn into the region r < r ; the degree to which they are drawn into this 
region is dictated by the value of 1/A at the outer minimum at criticality. The smaller that 
minimum value of 1/A, the more contained the matter fields are. Since they have nontrivial 
fields outside the horizon the critical solutions for large b are examples of extremal black 
holes with non-Abelian hair. 

The qualitative picture differs somewhat when b ^ 400, in that the 1/A always has only 
a single minimum. One may think of the sequence of monopole solutions here as large-6 
solutions in which the inner minimum of 1/A drops out sufficiently early that no double 
minimum solutions exist. To illustrate this, in Fig. |7] we show the approach to criticality for 
b = 1000. 

D. Behavior of the critical solutions 

The critical solutions themselves are not accessible through our numerical method. Nev- 
ertheless, we can obtain good approximations to these by following sequences of regular 
monopole solutions. Here we briefly summarize how some of the properties inferred from 
these numerical solutions vary with b. 

As we have already noted, the critical solutions for small b have horizons at the Reissner- 
Nordstrom radius r , while for large b the horizon occurs at a value r* < r (see Fig. [|). 
There is a discontinuity in passing from the low-6 to the high-6 regime: The critical solution 
does not continuously vary from one type to the other, but instead undergoes something 
like a first-order transition. For a value of b just below the transition, an inner minimum 
appears^] at a radius r < tq, while 1/A has a double zero at tq. As b increases, the inner 
minimum rapidly descends and at the transition 1/A has two degenerate minima, one at 

< ro and the other at r$. As b increases further the outer minimum moves upward, while 
the inner minimum at r* remains at zero. Eventually the outer minimum vanishes, and the 
only minimum in 1/A is associated with the horizon at r*. This radius r* asymptotes to 
zero as b —>■ oo. As b is increased, the critical black hole solution approaches that found in 
the literature in the infinite b limit ||.[5|]. 



2 This inner minimum in low-6 monopole solutions only appears for b very close to the transition 
point and when a is very near a cr . 
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FIG. 7. Monopole solutions for b = 1000 and various values of a. The progression from 
dot-dashed line, to dashed line, to dotted line, to solid line corresponds to a = 0.5, 0.9, 1.0, and 
1.05576. The first three panels depict the metric functions (a) 1/A(r), (b) B(r), (c) (AB) l l 2 {r). 
(d) Matter fields variables; u(r) begins at unity at the origin and asymptotes to zero as r — > oo. 
h{r) is zero at the origin and asymptotes to unity as r — > oo. 
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FIG. 8. Extremal horizon radius for various values of b. Analytic arguments indicate that 
— > as b — > oo. 
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FIG. 9. Mass of critical solutions as a function of b. Low-6 type solutions (b < 40) all have 
= Mq. The dots indicate the values for high-6 type solutions. 
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FIG. 10. Critical values for a — 1 as function of b. The inset shows in detail the transition 
region between low-6 and high-6 behavior, with the open circle indicating the apparent transition 
point at b = 40. Note that at this point a = 1.4187 < 1.5, the latter value indicated by the dotted 
line. The dashed line is an extension of the high-6 type critical solutions to smaller b. For the 
solutions represented by this dashed line, the fields are well-behaved inside the horizon, but are 
not asymptotically flat at spatial infinity. 



The mass of the critical solution can be inferred from the long-range behavior of 1/A. In 
Fig. |9] we plot this mass as a function of b. As expected, the mass of the low-6 solutions is 
just that of an extremal Reissner-Nordstrom black hole. However, in the high-6 regime the 
mass of the critical solution is less than the extremal Reissner-Nordstrom value M . The 
mass decreases with increasing b, and appears to have an asymptotic value of about 0.990. 
Note that there is no discontinuity in the mass in going from the low-6 to the high-6 regime. 

Finally, Fig. |10| shows function of b. The critical value a cr is always of order unity 

and is a monotonically decreasing function of b. This can be understood by recalling that 
the mass of the flat space monopole increases (by a factor of 1.8) as the ratio of the Higgs 
mass to the vector mass varies varies from zero to infinity. Stated differently, the value of 
v needed to achieve a given mass decreases with increasing Higgs mass, making it plausible 
that the critical v should decrease in a similar fashion. 

For large b there appears to be a power-law relationship between b and a CT — 1. That 
a cr — > 1 as b —>■ oo is consistent with observations made in past investigations [|pp. The 
literature also || indicates that the curve should asymptote to a cr « 3.94 as b — > 0. The kink 
in the data at b ~ 40 is real and reflects the transition between the small and large b regimes. 
Note that a CT < 1.5 at this transition point; this fact will be of importance later. The dashed 
curve in the inset of Fig. |T0| shows an extension of the high-6 type critical solutions into the 
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\ow-b regime. The corresponding solutions are well behaved inside the horizon, but are not 
asymptotically flat at spatial infinity; see the Appendix for more details. 



IV. ANALYTIC CONSTRAINTS FROM EXTREMAL HORIZONS 

We would now like to use analytic methods to gain some deeper understanding of our 
numerical results. We begin by focusing on the critical solutions with extremal horizons, 
but regular at the origin, that are the limits of the gravitating monopole solutions. In this 
section we will obtain a set of conditions at the extremal horizon that are necessary, although 
not sufficient, for the existence of such solutions. These turn out to allow two distinct types 
of limiting solutions. We discuss both of these, as well as the nearby nearby nonsingular 
solutions, in more detail in the next section. Our analysis of the behavior of the fields near 
the horizon complements and extends earlier treatments HJTl|. 

Because a simultaneous zero of both If A and its first derivative is a singular point of the 
differential Eqs. (|2.8|) and (|2.9|), the critical solutions are nonanalytic at the horizon radius 



r*. Ordinarily, physical considerations would constrain the allowable singularities at such 
a point. However, we are not actually requiring that the solution be physically acceptable, 
but only that it be the limiting point of a family of physically acceptable solutions. Keeping 
this in mind, it seems reasonable to impose the following set of requirements: 

1. We assume that the functions u, h, If A, and B are all finite and continuous at r*. 

2. Since r is a singular coordinate at the horizon, we should not assume that the deriva- 
tives of these functions with respect to r are continuous, or even finite, at r*. Instead, 
we require only that d(l/A)/dr vanishes at the horizon and that du/dr and dhfdr 
diverge less rapidly than A 1 / 2 does. 

3. We assume that the leading singularities of these quantities near the horizon can be 
approximated by (not necessarily integer) powers of \r — r*|, although we allow for the 
possibility that both the power and the coefficient of the leading term may be different 
on opposite sides of the horizon. 



Given these assumptions, Eqs. (|2.8| ) and Q2.9|) imply that the matter fields at the horizon, 
which we denote by u* and h*, must lie at one of the stationary points of U that were 
enumerated in Sec. II. Once these are specified, r* is determined by Eq. ( |2.11| ), which reduces 
to 

= 1- 87iGrlU{rA . (4.1) 

The monopole solutions without horizons were solutions to a boundary value problem 
with conditions imposed at r = and r = oo. The existence of the extremal horizon imposes 
three more conditions (on the values of u, h, and 1/A at r = rA and thus leads to a pair of 
boundary value problems, one for the interval < r < r* and one for r* < r < oo, that must 
be solved simultaneously. The interior problem has three boundary conditions at the origin 
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and three at the horizon, for a total of six, while the exterior has two at spatial infinity and 
three at the horizon, giving five in all. 

A standard approach to such problems is to look for a family of solutions obeying the 
conditions at one of the boundaries. If these can be shown to depend on N adjustable 
parameters, integration of these solutions will (assuming no singularities intervene) give an 
^-parameter family of solutions in the neighborhood of the other boundary. Generically, 
a necessary, although not sufficient, condition for a solution is that N be greater than or 
equal to the number of conditions imposed at the second boundary. Thus, at first thought 
it would seem that, given a suitable set of values for -u*, h*, and r*, we would need a two- 
parameter family of solutions as r — r* — > + to solve the exterior boundary value problem 
and a three-parameter family for r — r* — > 0_ to solve the interior problem. However, we do 
not expect to find extremal solutions for arbitrary values of parameters. Instead, for fixed 
b we must adjust v (or, equivalently, a) to its critical value. We may therefore view v as an 
extra adjustable parameter, and need only require the existence of two-parameter families 
of solutions on both sides of the horizon. The nonanalyticity at the horizon allows us to 
choose the parameters on the two sides independently. 

There are several caveats. First, the presence of the appropriate number of adjustable 
parameters at the horizon does not guarantee the existence of a solution. Global consider- 
ations that are beyond the scope of this analysis may make it impossible to satisfy all of 
the boundary conditions. Second, it may, and in some cases does, happen that the leading 
behavior near the horizon is fixed and that the adjustable parameters appear only in sub- 
dominant terms. Finally, if the values of u and h at the horizon are the same as at the origin 
(infinity), the interior (exterior) boundary value problem for the matter fields has a trivial 
solution, and so it is not necessary to have any adjustable parameters. 

We find it convenient to define 



/ u — u* 
\er*v(h - h*)/\/2 i 



(4.2) 



and to use the dimensionless position variable x = (r — r*)/r*. It is also useful to define the 
matrix 



ij ~ 2 dijjidijjj 



(4.3) 



The orthonormal eigenvectors and the eigenvalues of M. play an important role in the 
analysis; we denote these by 9 a and ji a (a = 1, 2), respectively. Finally, we define the 
two-component vector 



e 2 5 Q2JJ 



■ i 



2 drdi/ji 
as well as the ratio 



(4.4) 
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a = 



4nG 

e 2 r 2 



(4.5) 



where, as before, r is the horizon radius of the extremal Reissner-Nordstrom black hole. 

Equation (|2.11| ) and the equations obtained by substitution of Eq. fl2.10|) into Eqs. (p.8|) 
and (|2.9| ) form a set of three equations for the functions u, h, and 1/A. In the notation we 
have just defined, the first two of these can be compactly written as 



.4 

while the third becomes 



\r + Q) V + W = + $x + • • • 



.4 



2Fx - 2a 



i) t Mi) + -{i)') t i)' 



+ 



where 



-ae 2 rl 



d{r 2 U) 



dr 



(4.6) 



(4.7) 



(4.8) 



Here primes denote derivatives with respect to x, while the ellipses represent terms, of higher 
order in either x or the components of ip, that can be neglected here. 

We are assuming that the leading behavior of the various functions can be approximated 
by powers of x. It is then fairly easy to showQ that 



1 

A 



kx 2 + 



This reduces Eqs.([4.6|) and fl4.7| ) to 



and 



k = F -a 



-ip l Mip + kx(ip'Yip' 
x 



+ 



(4.9) 



(4.10) 



(4.11) 



We now turn to the behavior of ip as x — > 0. We see that ip must vanish at least as fast 
as |x| 1//2 , since otherwise the nonlinear term on the left-hand side of the equation dominates 
and there is no solution. Hence, we may write 



3 Assume that 1/A ~ \x\ a . If a > 2, Eq. ( ^4.7|) is dominated by the two terms not involving 1/A, 
and one finds that i[> ~ |x| 1//2 . This in turn implies that Eq. (4.6) is dominated by a single term, 
the first one on the right hand side, and hence has no solution. If instead a < 2, the two terms 
involving 1/^4 dominate Eq. ([4.7]), again implying that if) ~ | a?| 1 / 2 . Equation ( [4.6[ ) is now dominated 
by the three terms on the left hand side. Because the last two of these cancel, there is again no 
solution. 
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if, 



Vl/2\X\ 1/2 + 



T)iX + r^x* -\ + ?7 7l |x| 71 + r/ 72 |x| 72 + 
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(4.12) 



where the ellipsis within the square brackets represents higher-order analytic terms, the 
7j > 1/2 are noninteger powers to determined, and the final ellipsis represents smaller 
nonanalytic terms that are determined by the lower-order terms. 

When this expansion is substituted into Eq. fl4.10| ), the terms of order |x| 1//2 give the 
nonlinear equation 



k 



a 



V1/2V1/2 



771/2 = Mrji/2 



(4.13) 



where e is equal to plus or minus unity according to whether x is negative or positive (i.e., 
for the interior and exterior problems, respectively). A solution is possible only if 771/2 is 
proportional to one of the eigenvectors of A4. Let us denote this eigenvector by 9\\ and the 
orthogonal eigenvector by 9±, with analogous conventions for the eigenvalues. We thus have 



V1/2 = p9\\ 



(4.14) 



with p obeying 



kp 



6 a 5 

e—p 

4 4 F 



V\\P 



In addition, substitution of our expansion for ijj into Eq. (|4.10|) yields 



a 



1 — e—p 



F + ecr/i||p 2 . 



(4.15) 



(4.16) 



These equations can be used to find fcasa function of p 2 , and the result then substituted 



back into Eq. (|4.15 ). For the interior solution (e = 1), this leads to the result that 



Pint 



2(T/i|| 



- 1 ± 




or 



Pint 



(4.17) 



and 



k 



4F 



int 



(2 



(4.18) 



Since k must be positive, we require that F > 0. The requirement that p be real then implies 
that nonzero solutions for p exist only if /xy > —F/4; for the solution with the lower choice 
of sign, there is the additional requirement that /in not lie between and 3F/4. 

Although Eqs. ( |4.15 ) and ( 4.16|) allow a nontrivial solution for p 2 in the exterior region, 
other constraints, described below, require that 



Pcxt 







(4.19) 



and hence that 
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k ex t = F 



(4.20) 



Having determined 771/2, we can now turn to the remaining terms in Eq. ( f4.12| ). Extract- 
ing the coefficients of the integral powers of x in Eq. ( 4.10| ) yields a series of inhomogeneous 
linear equations that determine the r] n ; the first of these isQ 



<7 r 



(2k-M)r]i-ekp 2 - 2 (ASt/i ) 0n + Vi 



(4.21) 



In a similar fashion, the nonanalytic Odxp) terms give a linear equation for t/ 7j . How- 
ever, because the inhomogeneous term in Eq. ( 4.10 ) is analytic and cannot contribute, the 
resulting equation is homogeneous. A solution of this equation is possible only if 7/ 7j is 
proportional to one of the 9 a and if jj is a root of 








7(7 + 1) 



e'yap 



7(7 + 1) - 2 e 7^ 



- AU 



V1/2 oc 6\\ 
r/1/2 oc 9_l . 



(4.22) 



Note that both equations take the same form in the p 2 = case, where the two \i a are on the 
same footing. Because the equation is homogeneous, 7/ 7 ,. is determined only up to an overall 
multiplicative constant. In order that the solution near the horizon have two adjustable 
constants, there must be two independent nonanalytic terms with powers greater than 1/2. 
Thus, the pair of equations corresponding to the two /i a must, between them, have two roots 
greater than 1/2. 

For p = 0, the coefficient of 7 in Eq. ( 4.22 ) is positive, implying that at least one 7 must 



be negative for each \i a . In order to have two roots positive and greater than 1/2, we must 
require that both eigenvalues satisfy 



3F n 
jl a > — > 



V 



. 



(4.23) 



For p ^ 0, the form of Eq. Q4.22|) depends on whether we are considering the mode 
proportional to 9\\ or the one proportional to 9±. In either case, it is most convenient to 
proceed by using Eq. ( |4.15| ) to eliminate /i a . For e = — 1 the coefficient of 7 is positive, so 
at least one root is negative for each mode. Detailed examination of the equation for the 9\\ 
mode then shows that it has no roots greater than 1/2; since the equation for the 9± mode 
can have at most one such, we must set p = in the exterior region. For e = 1 (i.e., the 
interior region), the equation for the 9u mode has one and only root greater than 1/2 for 



The coefficient of r/i on the left hand side of this equation vanishes for certain choices of pa- 
rameters, leaving r]i undetermined. A similar phenomenon can also happen in the equations for 
the other rj n . These parameter choices correspond to points where the power 7^ of one of the 
nonanalytic terms goes through an integer value. 



20 



all values of /xy consistent with the reality of p. The equation for the other mode can be 
rewritten as 

For this to have a solution with 7 > 1/2, either /ly < fi± or /jy < 0; since we will see that at 
most one of the /i a can be negative, the second alternative implies the first. Combining this 
with the previous conditions on the fi a , we see that a p 7^ solution exists in the interior 
region if 

H ± > fjL\\ > -- , p^O. (4.25) 



Equations ( 4.23|) and ( [4.25 ), together with the condition F > 0, are the fundamental 



conditions that must be satisfied at the extremal horizon. To explore the various possibilities, 
we must apply these in turn to the stationary points of U that were enumerated in Sec. II. 
The last two cases can be immediately eliminated: Case 3 is excluded because both \i a < 0, 
thus ruling out the possibility of a p = exterior solution, while Case 4 has F < 0. This 
leaves Cases 1 and 2. The former, = and h* = 1, corresponds to an extremal horizon 
in the Coulomb region while the latter, u* = u and h* = h, gives a horizon in the monopole 
core. 5 We will study these in detail in the next section. 

V. BEHAVIOR NEAR COULOMB AND CORE REGION EXTREMAL 

HORIZONS 

In Sec. IV we showed that there are only two possibilities for the values of the fields at 
the extremal horizon. Although the local analysis that we used cannot tell whether these 
are consistent with the existence of a global solution, our numerical results show that both 
types of solution actually occur. In this section we examine these more closely. We also 
consider the nonsingular near-critical solutions in which 1/A has a related minimum. 

A. Coulomb region horizons 

The first possibility, u* = 0, /t* = 1, corresponds to a horizon in the Coulomb region 
outside the monopole core. Equation (|4.1|) implies that r* = r , the extremal Reissner- 
Nordstrom value, so a — F — 1. The matrix M. is diagonal, with eigenvalues 

A t 1 = M n = (evr ) 2 - 1 = - - 1 



Since it seems unlikely that the minimum energy solutions will have either u or h change sign, 
we assume that the fields are positive at the horizon. There can, however, be excited monopoles 
in which u(r) goes through a zero; see, e.g., Ref. ||. 
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fl2 = M 
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2ba 



(5.1) 



A solution outside the horizon is given by the extremal Reissner-Nordstrom metric, with 
u = and h = 1 for all r* > r . We assume that this is the only exterior solution satisfying 
the boundary conditions, and concentrate on the interior solution. 

There are several possibilities for the behavior of fields just inside the horizon: 
Type I: The singularity in the matter fields is less singular than x 1 ^ 2 (i.e., p — 0). Equa- 
tion ( |4. 18| ) implies that k — 1, so near the horizon 



u = CJx| 7u + a u x + ■ 
h = l-C h 
1 2 

A= X+ -- 



x\ lh — ahX + 



(5.2) 



where the ellipses represent terms that are determined by the terms shown explicitly. Here a u 
and ah denote constants whose values are fixed by Eq. ( f4.21| ). C u and are not determined 
locally, but must instead be adjusted so the boundary conditions at the origin are satisfied. 
The exponents ^ u and 7^ are solutions of Eq. ( f4.22p and, depending on the values of a and 
b, may or may not be greater than unity. The requirement that these exponents both be 
greater than 1/2 leads to Eq. fl4.23|) , which implies that 



a > Max 



3 

8b 



(5.3) 



Type II: Matter fields with |a;| 1//2 singularities. There are two possibilities here, depending 
on which field has the singularity. In the first (Type Ha), 



u 

h 
1 



1 1/2 
C h x 



p\x\~'~ + C u |x| 7u + a u x + 



a h x 



J- = kx 2 + 
A 



(5.4) 



while in the second (Type lib), 



u 



C u \x\^ + 



h = l- p^fa\x\ 1/2 - C h \x\ lh - a h x + 

\ = kx 2 H 

A 



(5.5) 



Here 



(l-p 2 /2)" 



(5.6) 



and 



P 



2//|| 



4//|| - 1 ± Jl + 4//|| 



(5.7) 
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with equal to the lesser of the two eigenvalues in Eq. ( |5.1| ). If 2 < a < 7/2, the upper 
sign must be used in Eq. ( |5.7| ). For Type Ha solutions Eq. ( |4.25| ) requires that either 

3 2 1 

- < a < and 6 < - (5.8) 

2 1-46 4 v ' 



or 



3 1 
a > - and 6 > - . (5.9) 
2 4 v ; 

The corresponding requirement for Type lib is that 

2 1 

a > > 2 and b < - . (5.10) 

1-46 4 v ; 

At the boundary between the Ha and lib regimes, a = 2/(1 — 46), both 7^ and 7^ are equal 
to 1/2 and the two types of solutions merge. Note that none of these solutions is possible 
if a < 3/2; this point will be important when we compare this analysis to the numerical 
results. 

The most important difference between the Type I and the Type II solutions is in the 
behavior of the quantity yAB. If p ^ 0, one of the matter fields varies as |x| 1//2 near the 
horizon. As a result, the integral in the exponent in Eq. ( |2.12| ) diverges for any r < r , 
and a/ AB becomes a step function centered at the horizon. Horowitz and Ross [O have 
shown that a particle in geodesic motion across a horizon feels a tidal force proportional to 
the logarithmic derivative of V AB, and have used the term "naked black holes" to describe 
certain near-extremal black holes solutions for which this quantity is large. Not only do our 
Coulomb region extremal solutions fit in this category, but so do the nearby nonsingular 
solutions. 



B. Horizons in the monopole core 

The second possibility allowed by the analysis of Sec. IV is that w* = w(r*) and 
= h(r*), with u and h given by Eqs. ( |2.6| ) and Q2.7p . This corresponds to a horizon 
in the monopole core region. Unlike the case of a Coulomb region horizon, there must be a 
nontrivial exterior solution, thus giving an extremal black hole with Higgs and gauge boson 
hair. Because our numerical results show solutions of this type only for relatively large val- 
ues of 6, we will find it convenient to use large-6 expansions to simplify some of the algebra 
below. 

Substituting the expressions for it and h into the potential U gives 



8nGr 2 U = — -(evr) 



-1 + 2b(evrf - b(evrf . (5.11) 



Equation (4J.) requires that this be equal to unity at the horizon. By making use of the 



fact that (evr*) 2 = a/2a, this constraint can be rewritten as 
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The roots of this equation determine a in terms of 6 and a. For large 6, the 6-dependence is 
negligible, and 



a = 2a 



it vi 



a 



-i 



0(1/6) . 



(5.13) 



In contrast with the previous case, the matrix M. is not diagonal. Using the equations 
obeyed by u and h, it can be written as 



M 



2u 2 2y/2evr* uh 
2y/2evr, iih 4be 2 v 2 rlh 2 
( 2(1 -y 2 ) 



2v / 2y % /T^ 



V2 V / 2WT^F 4by 2 - 4(1 - y 2 
where in the second line we have defined 




0(1/6) 



(5.14) 



y 



tew*: 



2a 



(5.15) 



The eigenvalues of A4 are 



tlx = + 0(1) 

a 

^ 2 = 2- - + 0(1/6) . 
a 



(5.16) 



The exterior solution must be of the p = type, implying that these eigenvalues must 
obey Eq. ( 4.23 ). This requires that they both be greater than 3F/4, where 



6-1 

2 



a^b 
~4a 



a 



a 



+ 0(1/6) 



(5.17) 



In the large 6 limit, y^i clearly satisfies this condition, while for \xi we obtain the constraint 

(5.18) 



3 1 
a > —a 2 + -a + 0(1/6) . 
32 2 v 1 ' 



We also obtain the requirement 6 > 1, since otherwise M has a negative eigenvalue. Since 
a and a are also related by Eq. (|5.13|) , we can obtain a constraint that depends only on a. 
To satisfy both Eqs. ( |5.13|) and ( |5.18| ), we must take the upper sign in Eq. fl5.13Q . We can 
then combine the two conditions to obtain 



6a 2 - —a -4 + 0(1/6) > 
15 



(5.19) 
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which implies that a > 1.06 + 0(1/6). Note in particular that a = 1 (i.e., r* = ro) is not 
possible for any value of b. 

Once this condition is satisfied, there are no further constraints imposed by local analysis 
at the horizon. In particular, the interior solution can be either of the p = or the p ^ 
type. Because M. is not diagonal, as it was in the Coulomb case, the same irrational powers 
appear in both u and h. Thus, if p = 0, the fields near the horizon behave as 

u = u(r*) + g"Ci|x| 71 + g^C^M 72 + ol u x H 

h = h(n) + gfCi|x| 71 + q h 2 C 2 \x\^ + a h x + --- 

- = Fx 2 + • • • (5.20) 



Here 71 and 72 are solutions of Eq. fl4.22p corresponding to the two \i a \ as before, their values 



may be either greater or less than unity, depending on the values of a and b. The constants 
q% and q% are determined by the eigenvectors 9 a . The values of a u and oth are fixed and are 
the same on both sides of the horizon, while C\ and C 2 can be varied independently inside 
and outside the horizon so that the boundary conditions can be satisfied both at the origin 
and at spatial infinity. 

If instead p 7^ 0, the fields outside the horizon are as in Eq. ( |5.20| ), but inside the horizon 



= u(n) +pq%\x\ 1/2 + g^CiM 71 + q 2 C 2 \xp 2 + a u x + ■ ■ • 
= h(n) + pq 2 l \x\ 1/2 + gfCxIx) 71 + q\C 2 \x\^ +a h x + ■ ■ ■ 

Ut^ j . 



u 

h = h(r 

- = kx 2 + ■ ■ ■ (5.21) 

We have assumed that /i 2 < if the opposite were the case, the \x\ 1 / 2 terms would involve 
q\ and q\. Here k and p are given by Eqs. ( 4.18 ) and ( 4.17|) , with there being two acceptable 



solutions. As before, the values of a u and ah are determined, although they are no longer 
the same as in the external solution. 



VI. COMPARISON WITH NUMERICAL RESULTS 

Let us now see how well these analytic results agree with our numerical results. We 
begin with the low-6 regime, where the nonsingular monopole solutions tend toward a critical 
solution with a Coulomb region horizon at the Reissner-Nordstrom value r = r . Throughout 
this regime, the critical solutions we find appear to be the more singular Type II solutions 
rather than the Type I solutions of Eq. (|5.2| ). This is consistent with our finding (see 
Fig. |10D that the values of a cr are always less than the lower bound for Type I solutions 
given in Eq. (|5T5|). 

Given the values for function of b, Eqs. ( p.8[ ) and ( |5.10| ) predict that there should 

be a transition at b ~ 0.1 between the Type Ha regime, where u ~ |x| 1//2 and the Type lib 
regime, where 1 — h has the |x| 1//2 singularity. At b = 0.1, the lowest b for which we have a 
sequence of stable configurations that approaches a critical black hole solution, we find that 
the powers of u and 1 — h are both close to 1/2, as predicted for this transition. It is curious 
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to note that this is close to the value of b below which (1/A) m in can be a double- valued 
function a. 

Examining the Coulomb-type black hole solutions in more detail, we find, for b & 20, 
that as r approaches r from below, the values of d 2 (A~ l )/dx 2 and u 2 /x tend toward the 
values of k and p predicted by Eqs. ( |5.6| ) and ( |5.7| ) with the upper choice of sign. 

We observed in Sec. IIIB that the matter fields displayed exponential falloffs when written 
as functions of the proper length l(r) for low-6 type solutions. To understand this, we note 
that the field equation for u, Eq. (|2.8|), can be written as 



dl 2 + 



d(A- 1 / 2 ) 8%GrK 



+ 



dr J~A 



du u(u 2 - 1) 2 2 2 

— = 7, V e v uh L . (6.1) 

dl r z 



Our analytic results indicate that the quantity in square brackets is of order unity as x — > 
for the interior portion of the extremal solution. This suggests that this quantity should be 
a relatively slowly varying function for large I and hence that the falloff of u should behave 
as a decaying exponential with a slowly- varying decay constant. As we saw in Fig. |3|, this is 
indeed the case both for the critical and the near-critical solutions. A similar analysis can 
be applied to 1 — h. 

From this point of view, the vanishing of u and 1 — h at the extremal horizon is a 
simple consequence of the fact that l{ro) = oo. This analysis also tells us how these fields 
should behave as the critical solution is approached. We concentrate here on u, and restrict 
ourselves to the range of b where u ~ \x\ l l 2 in the critical solution. Near f, the minimum of 
1/A, we can writef] 

— ps c\(r - f) 2 + A 2 . (6.2) 

Substituting this into Eq. (|3.3| ), we find that 

I = 1(f) ~q 1 ln(l/A) . (6.3) 

Writing u as a decreasing exponential, u ~ u exp(— c 2 Z), we then have 

u(f) ~ u (l/AY min (6.4) 

with q = c 2 /2ci. This analysis does not give predictions for C\ and c 2 . However, we can get an 
indication of what to expect for q by recalling that for the critical (A = 0) case u ~ (l/A) 1 / 4 
as the horizon is approached from the inside. If the above expressions were exact and c\ and 
c 2 were truly constant and independent of A, this would imply that q = 1/4. Turning to 



6 Note that c 2 cannot be the same as the coefficient k that appears in the critical solution, since the 
second derivative of 1/A has a discontinuity at the horizon of the critical solution, but is continuous 
when A ^ 0. 
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our numerical results, we do indeed see a power law behavior for u, with q ranging between 
0.23 and 0.28. 

A similar analysis can be applied to the jump in V 'AB near the minimum of 1/A. For the 
parameter ranges where we find extremal solutions with p ^ 0, this predicts that vAB r =o 
varies as (l/A)^/ 4 assuming a normalization of a/ AB r=00 = 1. Not only does this account 
for the power law behavior noted in Sec. IIIB, but the predictions for the power are borne 
out by the data within the numerical accuracy, where the power p 2 /4 ranges from about 0.7 
to unity. 

Near b = 20, a cr equals 3/2, the lower bound allowed by the analysis of the previous 
section for Coulomb-type solution. Nevertheless, the qualitative nature of the approach to 
criticality does not change as b increases past this value and as a cr becomes smaller than 3/2. 
Indeed, the transition to core-type horizons does not appear to occur until about b = 40, 
where a CT = 1.4187. It does not seem as if this discrepancy can be attributed simply to 
numerical errors. We find that changing numerical parameters or modifying details of the 
algorithm gives a variation in the value of a cr at the transition that is at most of order 0.01. 
For more details of the numerical aspect of this issue, see the Appendix. 

The resolution of this puzzle seems to lie beyond the precision of our numerical simula- 
tions. Several possibilities suggest themselves. First, there may be a sharp change in the 
approach to criticality that is only seen when one reaches solutions for which (l/74) m i n is 
less than the limits set by our spatial step size, in such a way that our extrapolation to the 
critical limit is incorrect. It could be that the decrease in the minimum of 1/A suddenly 
slows, and that it only reaches zero at a cr > 3/2, or it could be that a new family of solutions 
intervenes. Alternatively, it may be that a cr is correctly extrapolated, but that the critical 
solution violates one or more of the assumptions that we enumerated at the beginning of 
Sec. IV. For example, the critical solution might have discontinuities in the matter fields at 
the horizon; because this would occur at a zero in 1/A, it would not necessarily lead to a 
divergent energy. 

We turn now to the high-6 regime, where the critical solution has an extremal horizon 
in the core region, with = u(r*), = h{r*). The numerical solutions all have p = 
for both the internal and external solutions, with the matter fields varying less rapidly than 
| re | 1//2 near the horizon. Moreover, for the entire region where we observe this type of black 
hole solution, the leading behavior of u — and h — at the horizon appears to be linear 
in x, with the nonanalytic terms being subdominant. The relations between a, u*, h* and 
r* are all in agreement with our analytic results. 

Examining the solutions in more detail, we find that in most of the interior portion of 
the critical and near-critical solutions both u(r) and h(r) closely track u{r) and h(r). This is 
easily understood by recalling the form of the matter action, Eq. ( |2.5| ). Once 1/A has become 
small, the gradient terms in the action are much less important than the position-dependent 
potential U . As a result, the fields that minimize the action should point-by-point be close 
to the minimum of U, namely u{r) and h(r). 

It may be at first surprising that a second minimum in 1/A exists at r pa tq outside the 
true horizon. To explain this, consider the behavior of the metric outside the monopole core. 



27 



If the matter fields u and h are negligibly small for r > R COIC , then in this Coulomb region 
the mass function for a configuration with magnetic charge 4n/e can be approximated by 



m{r) = M 




(6.5) 



This gives a local maximum of m(r)/r, corresponding to a local minimum of 1/A, at 



with the value of 1/A at this minimum being 1 — (M/Mq) 2 . This minimum will occur, 
regardless of whether or not there is a horizon at a smaller value of r, as long as r min > -R cor c- 
It is absent in our solutions for very large b because for those solutions a, and thus v, are 
small enough that the core extends beyond r . 



We have used both analytic and numerical methods to study self-gravitating Yang-Mills- 
Higgs magnetic monopoles for a range of parameters, with an emphasis on their approach to 
the black hole threshold as the Higgs expectation value tends toward its critical value. We 
find two quite distinct behaviors. In the low-6 regime, with weak Higgs self-coupling, the 
critical solution is identical to the extremal Reissner-Nordstrom solution outside the horizon. 
Inside the horizon, the matter fields are smooth and nonsingular, as is 1/A. However, because 
of the step function behavior of V AB at the horizon, B = g tt vanishes identically in the 
interior. The associated singularity at the horizon means that an observer falling freely into 
a critical or near-critical monopole will experience the strong tidal forces characteristic of a 
Horowitz-Ross naked black hole. 

One might speculate that these singularities were inextricably associated with the tran- 
sition from a nonsingular monopole to a black hole spacetime and were independent of the 
details of the Higgs potential. However, our solutions in the large-6 regime show that this is 
not the case. For these, the critical solution has nontrivial matter fields outside the horizon, 
and thus is an extremal black hole with hair. Although the metric functions are still non- 
analytic at the horizon, their singularities are much weaker and do not lead to large tidal 
forces. 

The transition between these two regimes is itself quite interesting, and is in many 
ways reminiscent of a first-order transition. For a range of values of b on both sides of 
the transition point, 1/A for the critical solution has two distinct minima. For values of b 
below the transition, 1/A is positive at the inner minimum and has a double zero at the 
outer minimum at r = r Q . The inner minimum moves downward with increasing b until, at 
the transition point, there are two distinct but degenerate minima, with 1/A vanishing at 
both. As b is increased further, the outer minimum moves upward while the inner minimum 
remains a zero of 1/A. While this implies a discontinuity in the horizon area, the mass of 
the black hole is continuous. In terms of black hole thermodynamics, this corresponds to 




(6.6) 



VII. CONCLUDING REMARKS 
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a discontinuity in the entropy but (because these are zero temperature black holes) a free 
energy that is a continuous function of b. 

Our focus in this paper has been on the behavior of the monopole solutions as the Higgs 
expectation value is increased to its critical value. It seems natural to ask how these solu- 
tions evolve if v is increased still further. Actually, there is a bifurcation at this point. If one 
requires that the solutions remain well-behaved at spatial infinity, the solutions with v > v CI 
are black hole solutions with singularities at r = 0. In the low- 6 regime, these are simply 
nonextremal Reissner-Nordstrom black holes, of varying mass, with no dependence on either 
a or b. (There are also magnetically-charged black holes with hair in this regime; however, 
they are not continuously connected to our critical solutions.) In the high-6 regime, the 
counting of boundary conditions suggests that, for given values of a and b, the continuation 
of the external solution is a one-parameter family of black holes with hair. Alternatively, 
one can require that the fields remain nonsingular at r = and continue the interior solu- 
tion. Because there are three boundary conditions at the origin, compared to the two at 
spatial infinity, this gives solutions with nonsingular interiors bounded by horizons that are 
uniquely determined by specifying a and b. The behavior of the metric near the horizon is 
similar to that found when de Sitter spacetime is written in static coordinates, suggesting 
a cosmo logical interpretation for these solutions. We suspect that, possibly with a suitably 
modified Higgs potential, these may give rise to topological inflation [|13| , p!4| . 

Finally, having studied these self-gravitating monopoles and obtained an understanding 
of their approach to criticality, we can use these solutions as tools for investigating the 
transition from a nonsingular spacetime to one containing blackhole horizons. We will 
describe this elsewhere. 
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APPENDIX A: NUMERICAL DISCUSSION OF TRANSITION REGIME 



In Section VI we noted that a discrepancy appears between the analytic argument that 
low-6 type critical solutions can only exist when a cr > 1.5 and the numerical observation 
that this type of critical solution persists when a CT < 1.5. We discuss here the details of our 
investigation that lead us to believe that this discrepancy is not a numerical artifact. 

Recall that the transition between the high and low-6 type critical solutions is reminiscent 
of a first-order phase transition. The metric function 1/A(r) for a regular monopole solution 
near criticality with b near the transition value possesses two minima, one near the extremal 
Reissner-Nordstrom radius, t~o, the other in the core of the monopole, r < r$. As a is 
increased (for a fixed 6), one of these minima becomes a horizon, and a critical solution is 
attained. If 1/A at the inner minimum becomes zero first, then we have a high-6 type critical 
solution; when 1/A for the outer minimum touches down first, we have low-6 type solution. 
Let us examine each of these minima separately for regular solutions near criticality in this 
transition region. 

1. The core minimum 

The inner minimum for the types of solutions under consideration varies rapidly with 
both a and b, so one can identify very precisely when such a minimum becomes a horizon. 
Moreover, the field variables and their derivatives are well-behaved near this radius, and the 
numerical solutions here are robust and insensitive to variations in the numerical algorithm. 
Different numerical implementations of the field equations leads to differences in a cr of a few 
parts in 10 4 for a high-6 type critical solution. 

The numerical and analytic results could be reconciled if (1) the high-6 type solutions 
persisted down to lower b into a region where a cr > 1.5 and (2) the identification of the 
outer minimum as a zero of 1/A for the critical solutions in this region was a result of 
numerical error. We can address the first issue by using a trick. In Sec. IV we learned 
that critical solutions solve separate interior and exterior boundary value problems. Hence, 
we can determine function of b for the high-6 type solutions by just solving the 

field equations from the origin to the horizon, ignoring the exterior region from the horizon 
out to spatial infinity. By using an algorithm similar to that described in Sec. Ill, we find 
the family of critical solutions represented by the dashed line in the inset of Figure [1(| 
These solutions have high-6 type horizons but are not asymptotically flat at spatial infinity; 
however, if there were an asymptotically flat critical solution with the same value of b, it 
would have the same a cr as these solutions. This dashed line crosses a cr = 1.5 when b = 25.6. 
Thus to be consistent with the analytic results, the transition value between low and high-6 
critical solutions must occur for b < 25.6 (rather that at the value of b = 40 indicated by 
our numerical results). Again, this result appears to be numerically robust to a few parts in 
10 4 . 
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2. The Coulomb minimum 



According to the previous discussion, the ability to reconcile numerical and analytic 
results now depends on the sensitivity of the outer minimum, the one near r , to numerical 
error. From the inset in Figure [TJj, we see that in order for b tr < 25.6, our values for a cr for 
the low-6 type critical solutions must be in error (too low) by approximately 0.04. A priori, 
such an error is certainly plausible. At this minimum, the solutions near criticality vary 
sluggishly with a and b. The spatial derivatives of field variables are large at this radius, 
and indeed are expected to be singular if the outer minimum becomes a horizon. It is here 
that we are most susceptible to numerical errors. 

As we discussed in Sec. Ill, discretized versions of Eqs. ( p.8|) -( ^TT|) are implemented 
in our numerical algorithm to identify static monopole solutions. In order to probe the 
robustness of the algorithm and the confidence with which we can treat these solutions, 
different implementations of the discretization were examined and compared. In particular, 



we focus on Eqs. 
equations. 



and ( |2.9|) . First, we used the straightforward expansion of these 
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In a separate treatment, we used Eq. (|2.11|) to eliminate d(l/A)/dr. 
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Though these two systems of equations are equivalent in the continuum limit, their dis- 
cretized forms which appear in the numerical implementation suffer a (9[(Ar) 2 ] difference. 

Figure [11] shows the (1/A) min as a function of a for solutions near criticality using the 
two different equation sets Eq. ( |Al| ) and ( |A2|) at b = 16 and b = 35. The curves in Fig. |11 
do not extend into the region where (l/v4) min < 1.0 x 10~ 6 since this is 0[(Ar) 2 ] and our 
solutions break down. (For these calculations we used Ar = 5 x 10~ 4 .) Indeed we find 
no stable solutions at all in this regime; however, we expect this effect is an artifact of 
the discretization. What happens for this very small (1/A) min is beyond the ability of our 
resources to probe. 

Nevertheless, the extrapolated difference in a cr between the two numerical implemen- 
tations were on the order 0.001. Other similar variations in the numerical algorithm were 
tried, and analogous results are found. Moreover, varying the spatial discretization size also 
yielded differences in a cr on the order 0.001. Thus, numerical error does not seem to provide 
the Aa cr ~ 0.04 needed to account for the discrepancy between the numerical observations 
and the analytic results. To illustrate the scope of the discrepancy needed, the difference be- 
tween the two curves needs to be approximately eight times the total range of the horizontal 



axis shown in Fig. 11 



Finally, there is no apparent qualitative change between solutions whose critical a cr is 
greater than or less than 1.5, as seen in Fig. [TT]. The approach of (1/A) min to zero is a much 
slower function of a as b increases, but there is no sharp change at a cr = 1.5. 
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FIG. 11. Comparison of (l/^4) m in versus a using two different discretizations at b = 16 (top 
diagram) and b = 35 (bottom diagram). The dotted curve employs the system Eqs. (Al) while the 
solid curve employs the system Eqs. (A2). 
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